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ABSTRACT 

Feedback by active galactic nuclei (AGN) is frequently invoked to explain the cut-off 
of the galaxy luminosity function at the bright end and the absence of cooling flows in 
galaxy clusters. Meanwhile, there are recent observations of shock fronts around radio- 
loud AGN. Using realistic 3D simulations of jets in a galaxy cluster, we address the 
question what fraction of the energy of active galactic nuclei is dissipated in shocks. 
We find that weak shocks that encompass the AGN have Mach numbers of 1.1-1.2 and 
dissipate at least 2% of the mechanical luminosity of the AGN. In a realistic cluster 
medium, even a continuous jet can lead to multiple shock structures, which may lead 
to an overestimate of the AGN duty cycles inferred from the spatial distribution of 
waves. 
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1 INTRODUCTION 

There is increasing evidence from deep X-ray observa- 
tions for shock waves from radio-loud active galactic 
nuclei (AGN) in the ce ntral cores of galaxy clusters (e.g. 



iMcNamara et aO (|2005|)) . Exam p les are the (weak) shocks 
in M87 ( Formanet aU |2005| . ISimionescu et all l200'iD 
Hydr a A pulsen et all |2006| ) and Perseus (IFabian et al.l 
20061). Using d eeper Chandra data of the Virgo cluster, 
Forman et al.l (|2006l ) confirmed the presence of a weak 
shock at 14 kpc and determined its Mach number as 
M ~ 1.2. This shock was c onfirmed spectros co picall y 
using XMM-New t on d ata by ISimionescu et aL I (|2007l ). 
ISanders fc Fabian ('2006') report isothermal shocks in Abell 
2199 and 2A 0335-f 096 with Mach numbers of ~ 1.5. Shock 
waves have also been detected in the per iphery of the 
low-p ower isolated radio-galaxy NGC 3801 (jCroston et al.l 
I2OO7I ) and a strong shock has been associate d with the 
expa nding FR I radio galaxy in Centaurus A (|Kraft et al.l 
I2OO3I ). 



The detection of shocks around FR I sources has led 
to the suggestion that most, if not all, radio galaxies go 
through a phase that is associated wit h shock heating by a 
supersonically expanding radio source. ICroston" et al.l (l2007l ) 
have estimated that the energy stored in the shocked shell 
is equivalent to the thermal energy within ^11 kpc of the 
galaxy centre and a factor 25 larger than the inferred p dV 
work required to inflate the radio lobes. This suggests that in 
the early phases of radio-source evolution, the energy trans- 



fer from the AGN to its environment is dominated by shock 

hea ting. 

iNulsen et all l|2006D made the following estimate for the 
heat input into the ICM by shock waves: The heat per unit 
mass generated by a shock is given by 



AQ ^TAS = EAlnp/p'' 



(1) 



where E is the specific thermal energy, p the pressure 
and p the density of the gas. Thus, the fraction of the 
thermal energy that is dissipated, AQ/E, is given by the 
jump of In p/p'' in the shock. Three weak shocks are visible 
in the X-ray image of M87. For the innermost shock at ~ 
3.7 kpc a Mach number of 1.4 has been inferred, which 
implies a heat input of AQ/E « 0.022 and a shock age of 
2.4 X 10^ yrs. Obviously, the heat input of this shock is tiny. 
However, two more shocks have been identified at larger 
radii that require several times more energy. Thus, a shock 
of comparable strength to the 3.7 kpc shock may well occur 
every ~ 2.5 x lO'' y. The cooling time of the gas at 3.7 kpc 
is ~ 2.5 X 10* yrs, so that there is time for ~ 100 such 
shocks during the cooling time. Therefore, the combined 
heat input from ~ 100 of these shocks is more than enough 
to make up for radiative losses from the gas. 

In this Letter, we investigate what fraction of the jet 
energy is dissipated in shocks around the supersonically ex- 
panding radio source. Using a hydrodynamical simulation 
of jets in a realistic cluster set-up and a shock finding algo- 
rithm, we quantify the properties of the shock and the effect 
on the intracluster medium. 
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2 METHOD 

The initial conditions of our simulation are based on a rerun 
of the S2 cluster from Springcl ct al. (2001), whose proper- 
ties are sufficiently close to a typical, massive, X-ray bright 
cluster with a mass of A/ ~ 7 x IO^^'M© and a central tem- 
perature of 6 keV. The cluster appears as a classical, re- 
laxed cooling flow cluster in X-rays. Its density rises steeply 
in the centre, and the profile is very simila r to the den- 
sity profiles reported in Vikhlinin et al.l (|2006| ). The setup is 
the same as used in iHeinz et all (|2006| ). The output of the 
GADGET SPH simulation serves as the i nitial conditions for 
our simulation. We use the FLASH code l|Frvxell et al.ll20o"ol ) 
which is a modular block-structured adaptive mesh refine- 
ment code, parallelised using the Message Passing Interface. 
It solves the Riemann problem on a Cartesian grid using the 
Piecewise-Parabolic Method. Our simulation includes 7x 10^ 
dark matter particles. For the relatively short physical time 
of the jet simulation (25 Myrs), radiative cooling and star 
formation are neglected, though they were included in the 
constitutive SPH simulation. 

The computational domain is a 2.8 Mpc"^ box around 
the cluster's centre of mass. The maximum resolution at 
the grid centre corresponds to a cell size of 174 pc, implying 
11 levels of refinement. The simulations presented in this 
letter were performed assuming an adiabatic equation of 
state with a uniform adiabatic index of 7 = 5/3. 

The jet is injected through a nozzle placed at the cen- 
tre of the gravitational potential, coincident with the gas 
density peak of the central elliptical galaxy. The nozzle is 
modeled as two circular back-to-back inflow boundaries 2 
kpc or 12 resolution elements in diameter. The nozzle faces 
obey inflow boundary conditions fixed by the jet's mass-, 
momentum-, and energy fluxes. This treatment avoids the 
entrainment of cluster gas into the jet which is unavoidable 
in simpler schemes where the jet is approximated by inject- 
ing mass, momentum, and energy into a flnite volume of the 
cluster that contains thermal gas and is part of the active 
computational grid. We were thus able to separate cleanly 
jet fluid and cluster fluid in order to study the heat input 
into the ICM only. The jet is centred on the gravitational 
potential of the cluster and follows the (slow) motion of the 
cluster through the computational domain. 

The jet material is injected equally in opposite direc- 
tions with velocity Vjot = 3 x lO^cms"^ and an internal 
Mach number of 32. The jet power of the simulation pre- 
sented in this letter was chosen to be Wjct = 3 x 10*^ ergs s^^ , 
corresponding to a rather powerful source. Comparing this 
luminosity w i th th e sample of cavity systems studied by 
iBirzan et al.l (120041 '). the jet power in our simulation is at 
the extre me end of those ob served and well above that cited 
for M87 (|Allen et al.l ([2006')). However, even such powerful 
sources assume FR I morphologies in dense clusters such as 
this one (see, e.g., Perseus A). Even though Perseus A is 
inferred to ha ve a mean power output of ~ 10** erg s~^ 
(|Sanders fc Fa bian 2007), its peak luminosities are likely to 
be significantly larger. Hercules A has a power of 1.6 x 10*® 
erg s~^ as implied by the large-scale shocks found around 
it (|Nulsen et al. I I2OO5I : iMcNamara et al.ll2005l ) but still ex- 
hibits an FR I/II morphology. In our simulation, we chose 



such a high luminosity to ensure that the jet is able to push 
through the dense gas of the central galaxy. 

The shocks in our simulation are detected using a multi- 
dimensional shock detection module adopted from the sPPM 
codeQ based on pressure jumps across the shock. The basic 
algorithm evaluates the jump in pressure in the direction of 
compression (determined by looking at the velocity field). 
If the total velocity divergence is negative and the relative 
pressure jump across the compression front is larger than 
some chosen value {Ap/p ^ 0.25), then a zone is marked as 
shocked. Using the jumps across the shock in the 3 velocity 
components, we get the x-, y-, and 2-components of a unit 
vector pointing in the direction of the velocity jump, hence 
in the direction normal to the shock front. We now project 
the pre- and post-shock velocities onto the shock normal. 
The upstream and downstream pressure, the upstream ve- 
locity and density are then written out. We have tested this 
shock detection algorithm with one- and two-dimensional 
shock tube problems and found that the jumps in pressure 
and density are reproduced very well. The shock structures 
in a slice through the central regions of our computational 
domain are shown in Fig. [T] 



3 RESULTS 

We have simulated a jet that resembles a FR I source in 
a realistic cluster environment. Our simulation reproduces 
the shock structure in the inner 100 kpc of the cluster 
around an AGN with FR I morphology. One can see how 
the jet inflates bubbles that break off and start to rise 
through the cluster medium. This morphology resembles 
many of the low-power AGN that are observed to inflate 
bubbles at the centres of cooling flow clusters. 

Fig. [1] reveals two kinds of shocks: At the working 
surface of the jet, i.e. at the point where the jet impacts 
the ambient medium, there is a very strong shock wave. 
In the first million years after the start of the jet, this 
shock wave has Mach numbers with respect to the ICM of 
> 30. The normal of this shock surface is equal or close to 
the direction of the jet and does not encompass the entire 
jet region, contrary to the outer shock that is described 
below. Later on, as the ICM near the jet gets hotter, the 
Mach number of the strong inner shock decreases to close 
to ~ 8 — 10. When motions of the ambient medium cause 
the jet to break off or to change direction, this shock can 
detach and a new working surface forms. As the jet jitters 
and ambient material moves into the jet, multiple shock 
fronts develop in the region close to the AGN. This leads 
to multiple shock features, as one can see in the botton left 
panel of Fig. [T] In our simulation, we see that a continuous 
jet can lead to multiple shock fronts as those observed in 
M87 and other clusters. Hence, the exitence of multiple 
shock fronts does not necessarily imply an intermittency of 
the AGN. 



1 S.E.Anderson and P.R.Woodward, World Wide Web 
|http: / /www. lose. umn.edu/research/sppm Laboratory for Com- 
putational Science and Engineering, University of Minnesota 
(1995). 
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Figure 1. Slices tlirougli the cluster centre showing the shocks (left) and the gas density (right) at 5 Myrs and 15 Myr after the start of 
the jet. Shown is only the central part of the computational domain. The entire computational domain represents a volume of 2.8 Mpc'^. 



The second kind of shock is a weak and nearly spherical 
shock that travels from the point of the injection region 
outwards through the cluster. The Mach number of the 
outer shock remains at fairly constant values of around 1.1 
- 1.2 for the largest part of its propagation through the 
core of the cluster. The pressure jump across a shock of 
M = 1.2 is 1.55. The outer shock is a pressure wave that 
is driven by the additional pressure from the injected gas 
in the core of the cluster. This is different from the strong, 
inner shock that is driven by the ram pressure of the jet. 
After about 10 Myrs the outer shock becomes prolate in the 
direction of the jet. As the bubbles rise mainly in the jet 
direction, the pressure also increases preferentially in the di- 
rection of the jet. The Mach number is slightly higher in the 
direction of the jet than at the sides of the outer shock front. 

Next, we wish to compute the total energy thermalised 
in the shock front. One can write for the thermal energy flux 
generated at the shock: 

F = [ed - eu{pd/pur]vd, (2) 
where the subscripts d and u denote down- and upstream 



quantities, respectively, e is the thermal energy density, p 
gas density and v velocity. The second term inside the brack- 
ets subtracts the effect of adiabatic compression suffered at 
a shock. The total thermalised energy input per time (by 
shocks) divided by the mechanical luminosity of the jet is 
shown in Fig. [2] About 2% of the mechanical luminosity 
of the jet are converted lastingly into internal energy. The 
properties of the outer shock are found to be relatively in- 
sensitive to the mechanical luminosity of the jet. The inner, 
strong shocks are much more efficient at generating energy 
because their Mach number is much higher. However, their 
area is relatively small and thus they may be important for 
the interstellar medium of the host galaxy, but they are un- 
likely to have a significant effect on the thermal state of the 
ICM. The properties of the inner shock depend also quite 
sensitively on the exact jet parameters and to some degree 
on the numerical resolution of our grid. Meanwhile, the prop- 
erties of the outer shock are not sensitive to the numerical 
resolution and appear converged. We note that the timescale 
and geometry of the initial energy release is important for 
the computation of the energy deposited in the ICM or in- 
terstellar medium. If the same total energy were injected in 
the form of thermal energy in pressure equilibrium with the 
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surroundings, the amount of energy transferred to the ICM 
in the form of shocks would be very different. Even taking 
into account the uncertainties of our jet model, the approach 
presented here is much more realistic than schemes in which 
the energy is injected in pressure equilibrium. 

The total increase in internal energy per unit time of 
the ambient medium (i.e. excluding the jet material) within 
the outer shock divided by the total mechanical power of 
the jet is shown in Fig. [S] We see that for the first 15 Myrs 
of the AGN activity, a bit more than 30 % of the injected 
energy has been converted to thermal energy in the inner 
core. This is much more than what has been thermalised 
by the outer shock. The temperature increase in the cluster 
core is mainly caused by pdV work from the expanding 
bubbles. While the shocks thermalise only a few per cent 
of the jet energy, this also raises the entropy of the cluster. 
On the other hand, pdV work by the expanding bubbles 
is adiabatic until the gas motions induced by the rising 
bubb les are dissipated by viscous processes (|Nulsen et al.l 
I2OO6I ). Yet, a succession of shocks can be sufficient to offset 
the radiative cooling of the entire ICM, as only 4-5 shock 
fronts that are permanently present lead to a conversion of 
~ 10 % of the jet energy. 

Clearly, magnetic fields can al ter the dynamics o f the je t 
and the radio lobes. As shown in iRuszkowski et al.l (|2007h . 
large-scale cluster magnetic fields tend to drape themselves 
around the rising bubbles and can suppress the fragmenta- 
tion of the bubbles. The ensuing dynamics can be different 
from the dynamics modelled here, but the outer shock struc- 
ture is unlikely to be substantially different. In the presence 
of a physical viscosity with values close to the Spitzer value, 
the dissipation of the outer, weak shoc ks can be higher than 
what we find in ou r simulation (see also lBriiggen et aLll2005l : 
IRuszkowski et al.| [2004V 

We have presented a single simulation of a shock that 
has been produced by a FR I jet. Obviously, there are a lot 
of parameters that can be varied. However, these kind of 
simulations are very expensive; on NAS Columbia a single 
simulation took more than 12,000 CPU hours (on 64 proces- 
sors). Hence, we chose one exemplary case that reproduces 
properties that match observations such as those in Hercules 
or M87 (despite wide discrepancies in jet powers). In the long 
term, a proper parameter study ought to be conducted. 
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Figure 2. Plot of the energy generated at the outer shock surface 
per time divided by the mechanical luminosity of the jet. 



§ 0.4 




5 10 15 20 

time / Myr 



ACKNOWLEDGEMENTS 

The anonymous referee is thanked for a helpful report. MB 
and ER acknowledge the support by the DFG grant BR 
2026/3 within the Priority Programme "Witnesses of Cos- 
mic History" and the supercomputing grants NIC 2195 and 
2256 at the John-Neumann Institut at the Forschungszen- 
trum Jiilich. SH acknowledges support through NASA grant 
TM5-6007X. We also acknowledge a supercomputing grant 
on NAS Columbia. 

The results presented were produced using the FLASH 
code, a product of the DOE ASC/ Alliances-funded Center 
for Astrophysical Thermonuclear Flashes at the University 
of Chicago. 



Figure 3. Increase in internal energy of the ambient medium (i.e. 
excluding the jet material) within the outer shocks divided by the 
total mechanical energy injected so far. 
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